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Abstract Large spin systems as given by magnetic macromolecules or 
two-dimensional spin arrays rule out an exact diagonalization of the Hamil- 
tonian. Nevertheless, it is possible to derive upper and lower bounds of the 
minimal energies, i.e. the smallest energies for a given total spin S. 

The energy bounds are derived under additional assumptions on the 
topology of the coupling between the spins. The upper bound follows from 
"n-cyclicity" , which roughly means that the graph of interactions can be 
wrapped round a ring with n vertices. The lower bound improves earlier 
results and follows from "n- homogeneity" , i. e. from the assumption that the 
set of spins can be decomposed into n subsets where the interactions inside 
and between spins of different subsets fulfill certain homogeneity conditions. 
Many Heisenberg spin systems comply with both concepts such that both 
bounds are available. 

By investigating small systems which can be numerically diagonalized 
we find that the upper bounds are considerably closer to the true minimal 
energies than the lower ones. 

Key words spin system, XXZ model, Heisenberg model, energy bounds, 
partiteness 
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1 Introduction 

Rigorous results on spin systems like the Marshall-Peierls sign rule ^ and 
the famous theorems of Lieb, Schultz, and Mattis [2113] have sharpened our 
understanding of magnetic phenomena. In addition such results can serve 
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as a basis or source of inspiration for the development of approximate mod- 
els. For example, the inequalities of Lieb and Berezin relating spectral 
properties of quantum systems to those of their classical counterparts pro- 
vide a foundation for classical or semi-classical treatments of spin systems. 

In this article we will extend the body of rigorous results on Heiscnberg 
spin systems by generalizing the notion of "bi-partiteness" , which is funda- 
mental for the findings of Marshall, Lieb, Schultz, and Mattis PI^E). We 
will introduce two new concepts which rest on the topological properties of 
the interaction matrix connecting the spins of the systems. 

The first concept, n-cyclicity, uses the property of many spin systems 
that their "net" of interactions can be wrapped round an n-cycle. The tri- 
angular lattice may serve as an example, it can be mapped onto a triangle 
in an oriented manner as if one would wrap it round the triangle. In such 
cases an upper bound of the minimal energies -Bmin hi each subspace 7i(M) 
of total magnetic quantum number M can be derived for Heisenberg models 
and XXZ models in general. 

The second concept which leads to lower bounds rests on n-homogeneity, 
i. e. on the fact that the set of spins can be decomposed into n subsets 
of equal size where the interactions inside and between spins of different 
subsets fulfill certain conditions. 

Fortunately, many Heisenberg spin systems comply with both concepts 
such that both bounds are available. For all cases which were investigated 
it turns out that the upper bounds are rather close to the true minimal 
energies, whereas the lower bounds are not. Therefore, especially the upper 
bound can serve are a benchmark or guideline for approximate methods like 
DMRG or variational methods in order to rate the achieved quality. 

The resulting bounds improve earlier findings of Ref. {4,.^ji^ especially 
for frustrated spin systems. 

The article is organized as follows. In section |21 upper bounds and the 
concept of n-cyclicity will be discussed, in section |3 lower bounds and n- 
homogeneity will be introduced. Both sections start with subsections ex- 
plaining the idea followed by more mathematical subsections presenting the 
mathematical tools. At the end of each section the resulting bounds are 
given. Examples are provided in sectional A more technical calculation is 
carried out in the appendix. 

2 Upper bounds 

2.1 Idea 

It is obvious from the Ritz variational principle that an upper bound for the 
minimal energy can be provided if an appropriate trial state can be found for 
which the energy expectation value is known analytically and rather close 
to the exact ground state value. In the following we will construct such trial 
states for subspaces H{M), i. e. for total magnetic quantum number M. 
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Starting point is the magnon vacuum {M = Ns) which is mapped by means 
of suitable powers of the total ladder operator into the subspace TL{M). 
Since the total ladder operator commutes with the Heisenberg Hamilton 
operator this does not change the energy of that state. In a second step we 
assume a certain topological property of the spin array namely that it can 
be wrapped round an n-cycle and construct a generalized "Bloch operator" 
which is a unitary operator that adds appropriate phases to the components 
of the trial state. Utilizing the known action of the Bloch operator onto 
the Hamiltonian we can evaluate the energy expectation value analytically 
which results in the expression for the upper bound. 

In order to motivate our definitions in the next subsection we recall the 
definition of a bi-partite spin system in the case of constant coupling: It is 
required that the spin sites can be grouped into +sites and —sites such that 

only H — pairs are coupled but no ++ or pairs. We suggest the following 

generalization: Assume that complex phase factors e"^^ can be attached to 
spin sites j such that only constant phase differences |0j — (j)k \ occur between 
adjacent (coupled) spin sites. This is the requirement needed for the above- 
mentioned construction of the Bloch operator. The attachment of phase 
factors is no longer arbitrary if there are "loops" in the coupling scheme of 
the spin system, i. e. periodic sequences of adjacent spin sites. If only even 
loops exist we may choose the phase differences to be \<j)j — 4'k\ = and the 
system is bi-partite. However, in the case of odd loops it becomes necessary 
to "wrap" the loop around the complex unit circle and the resulting phase 
differences will be integer fractions of 2tt. We will make this more precise in 
the next subsection employing the language of graph theory. 



2.2 Definition of n-cyclicity 

In this section we consider systems with spin sites with spin s and con- 
stant anti-ferromagnetic coupling. Thus the complete information about the 
coupling scheme is encoded in some (undirected) graph 7 = (V, -T). The ver- 
tices of 7 are the spin sites, V = {1, . . . , N}, the set of edges of 7 consists of 
those pairs of sites which are coupled and will be denoted by F. We make 
the convention that («, j) S -T iff (j, i) G F and {i, i) ^ F. Hence the number 
of members of the set F, denoted by \F\, equals twice the number of bonds. 
Further, we will consider orientations on 7, denoted by 7^, i. e. we split F 
into disjoint subsets T = r+ U T", such that e F+ iff (j, z) G F~ . 

Then the Hamiltonian of XXZ-type can be written in the form 

('■j>er {t,j)er+ (jj>er- 

= A + g + g^, (2) 

where S > and denote the usual spin observables at site i with compo- 
nents s\'^\ ^ = 1,2,3, and sf = s[^^ ± isf \ Of course, only the splitting 
^ depends on the orientation, not the Hamiltonian itself. 



4 



Klaus Barwinkel, Heinz- Jiirgen Schmidt, Jiirgen Schnack 



In order to define a suitable concept of n-cyclicity we consider graph 
homomorphisms, i. e. maps between graphs, such that vertices are mapped 
onto vertices and the corresponding edges onto corresponding edges. Let Cn 
denote the cychc graph with n vertices which will be identified with the 
rt-th roots of unity 

£ = 0, . . . , n - 1 . (3) 

Further will denote the cyclic graph with anti-clockwise orientation. 

Any graph 7 (or the spin system itself) will be called n-cyclic or having 
the cyclicity n iff there exists a graph homomorphism 

: 7 — C„ . (4) 

In this case will induce an orientation on 7 in an obvious sense. 




Fig. 1 Tlie pentagon is 5-cyclic and also 3-cyclic (l.h.s.) whereas the tetrahedron 
is not 3-cyclic (r.h.s.), because if the four vertices of the tetrahedron are attached 
to the numbers 1,2,3 one number must repeat and occurs at adjacent vertices, 
which does not happen in the 3-cycle. 




It is only in certain cases that different cyclicities n and n' mean an 
essential distinctness. This is because for n > 4 any n-cyclic system is also 
[n — 2)-cyclic since three successive vertices and the corresponding edges 
can be mapped in a forward-backward-forward way, compare the l.h.s. of 
Fig. ^ which shows a homomorphism of a pentagon onto a triangle, as an 
example. Each 2r7i-ring and hence any 2m-cyclic system is n-cyclic for any 
positive integers m, n, since it is 2-cyclic and C2 can be homomorphically 
embedded into any n-cycle. 

Hence it makes only sense to distinguish between even-cyclic systems, 
which will be called 2-cyclic, and (2n-|-l)-cyclic system with maximal integer 
n. If a spin system is 2-cyclic in our sense it will be bi-partite in the sense of 
Rcfs. 12113) where, however, the theory also comprises cases with different 
coupling constants. 

We consider some more examples which illustrate the definition of cyclic- 
ity. A triangular plane lattice with suitable periodic boundary conditions is 
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3-cyclic, a square lattice or cubic lattice is 2-cyclic. The kagome lattice is 
3-cyclic but not 2-cyclic. 3-cyclicity is equivalent to 3-colorability. Hence the 
octahedron, the dodecahedron, the cuboctahedron, and the icosidodecahe- 
dron are 3-cyclic, cf. |7], but the tetrahedron is not, sec r.h.s. of Figure^ 

A natural basis for a matrix representation of H is provided by the 
product states | m ) = | mi . . . miq ), —s < rrii < s with 



s 

and 



(3)lm) =m,|m) , (5) 



s- \m.) = y/s{s + 1) - m,{mi ± 1) | mi, . . . ,m, ± 1, . . . , tun ) (6) 

for all i e {1, . . . , iV}. The state \Q) = | s, s, . . . , s) will be called the 
"magnon vacuum" . Further we define 

N 

a — '^Qi = Ns - M , ai ^ s - rrii for aU i G {1, . . . , A^} . (7) 



We also define a mapping h of product states into complex numbers which 
depends on the graph homomorphism by 

N 

h{ln)^l[h{^r . (8) 

1=1 

Then it is easily shown that if (m|tj|m') ^ then h{m) = e2'^*/"/i(m'). 

For any ^ = 0, . . . , iV — 1 we define a unitary "Bloch operator" (general- 
izing the corresponding definition for spin rings in Ref. 8 ) 

Ui : n{M) — > H{M) (9) 

by 

Ue\m) = h{mf\m) (10) 

and linear extension. Recall that ai — lixljn. Then the following relations 
hold: 

U\QUt = e-^^^g (11) 
UlHUe = A + cosaiiG + g^) - isinaeiG - g^) . (12) 

If -E'i„in(M) denotes the minimal energy eigenvalue within the sector T-l{M) 
and \lp) € 7i(M) is an arbitrary normalized state we have the obvious 
upper bound 

E^iniM) < . (13) 

The problem is to find a state \(p) such that {(f\H\(p) can be explicitly 
calculated and represents a close bound. To this end we map the magnon 
vacuum | i?) by (S^)'^ into TC{M), which remains an eigenstate of H with 
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the largest eigenvalue in the Heisenberg case ^ = 1, and change its phases 
according to the Bloch operator. More precisely, let 

\[2M)^CM{S-r\n) , (14) 

where Cm is the positive normalization factor, compare ()52|l . ensuring 
( ]7m I ) = 1 and define 

1^) = \ ^M,i) = Ui\Qm) , (15) 

Then we obtain 

As it must be, this result has the obvious value (5|-r|s^ for a = and remains 
unchanged under a <-> 2Ns — a. The proof of Eq. (|16|l is given in appendixIXl 
Now consider 

{ip\H\ip) = {f2M\ulHUe\f2M) (17) 
= {f2M\A + cosae{g + g^)\f2M) (18) 

= cosai{QM\^ + ig + g^)\f2M) (19) 



,^ cos ag . , ^ I . I ^ , 

+(1 ^){f2M\A\nM) 

= \r\s^ COS ai (20) 
_ cos a, 6\r\ / 2 _ 2sa{2Ns - a) 
^ S ' N \ '"^ 2Ns - 1 

In line (|18|l we used 1)12(1 and the fact that | i^A/ ) and G are real in the 
product basis of the |m) whence {f2M\{G — G'^)\^m) = 0. Eq. (^01) follows by 
Eq. 1161) and the observation that y + (^ + ^/^ ) is a Heisenberg Hamiltonian 
which has the eigenstate | SIm ) with eigenvalue |/^|s^. 

For spin rings and a — 1, \ip) is nothing else but the relative ground 
state. Generally for spin rings, | ip) has the same shift quantum number as 
the relative ground state [Hj- 



2.3 Analytical expression for the upper bound 

The best bound for E^in{M) is obtained if cos ae in (^0)1 is as low as possible, 
i. e. £ = ^ and cosa^ = —1 for even n and £ — ^iM for odd n. Therefore 
the upper bound is given by 

E^^{M) < c\r\s' + (1 - ^Ns' - ^^—^ j , (21) 
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where c = —1 in the case of even n and c = — cos — for odd n. Let (5 = 1 
and Ns be integer. Then the total ground state hes in the sector M = 0. In 
this case we obtain 

i?,ni„(0)<c|r|.^ + ^^^^, (22) 
which improves the upper Berczin-Licb bound Ef^ff^^^'^'^^s'^, see if c\r\ = 

TT^classical 
mill 



3 Lower bounds 

3.1 Idea 

For the lower bound to be derived in the following section the interaction 
matrix J = {Jfj,u) describing the coupling between spins at sites /i and v 
must have certain homogeneity properties. The matrix must be symmetric 
and must have constant row sum. This alone is sufficient to derive some 
lower bounds 'S', which can be improved using the topological structure of 
the interactions, as will be shown in the following. 

The derivation works by constructing another "averaged" Hamiltonian 
having an analytically diagonalizable interaction matrix J, which neverthe- 
less has only eigenvalues already present for the original interaction matrix 
J. Since, by construction, JJ > J, this relation also holds for the related 
Hamiltonians, and we arrive at a lower bound. Extending this idea the ob- 
tained lower bounds could be improved for particular systems. However, in 
this article we will confine ourselves to deriving bounds for general classes 
of systems. 

Using the topological structure of the interactions will further enable us 
to determine the degeneracy of some eigenvalues of J and therefore improve 
the calculations of Ref. 6 where this information was not exploited. 



3.2 Definition of n-homogeneity 

The Hamiltonian used in this section is different from that of the previous 
section and assumed to be of the form 

The matrix JJ of coupling constants J^j/ is assumed to be symmetric and 
having constant row sums j. The latter property can be viewed as a kind 
of gauge condition, since adding a diagonal matrix with vanishing trace to 
I does not change the Hamiltonian l^^ . see Ref. fTTij . 
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Being symmetrical, J has a complete set of (ordered) eigenvalues ji, . . . , jjv- 
One of them is the row sum j with 1 = "^(Ij Ij ■ • ■ > 1) as the correspond- 
ing eigenvector. Let JJ' denote the matrix J restricted to the subspace or- 
thogonal to 1, and jijiin the smallest eigenvalue of J', jmin may be to- fold 
degenerate. Further, we will denote the a-th normalized eigenvector of J by 

(Cla, . . . , CNa), 1- 6. 

^ ^ J ^u'^ua — ja^^a ; ^ ^ Cfia^fi^ — j i ^ — f ; ■ ■ ■ i ^ i (^^) 

where we also allow for the possibility to choose complex eigenvectors. Sums 
over a = 1, . . . , excluding aj will be denoted by where aj denotes 
the index (within the ordered set of all eigenvalues) of the eigenvalue j 
belonging to the eigenvector 1. 

For later use we will consider a transformation of the spin observables 
analogous to the transformation onto the eigenbasis of J and define 

T„EE^c^s^, andga = Tj;,-T„ , a = l,...,iV. (25) 

The inverse transformation then yields 

Sp = ^c^aT„, ^ = l,...,iV. (26) 

In particular, T^^ = S/Vlv. It then follows directly from the definitions 
that 

Ns{s + 1) = ^(s^)2 = ^ g„ ^ Is^ + ^ ' Q„ , (27) 



For a later use we also need a relation between Hamiltonians with different 
coupling matrices. Therefore, let H and H be two Hamiltonians of the form 
(|23|l with coupling matrices J and J, such that J > JJ. Then H > H. As usual 
the ordering ">" of operators is defined by the corresponding inequality 
for arbitrary expectation values. Since H depends linearly on J it suffices 
to show that JJ > implies H > 0. But this is obvious in view of (|28|l : 
H = Y^pj^Qp with if} > and Qfj > 0. 

Next we turn to the suitable definition of n-homogeneity. Let the set 
of spin sites {1, . . . , N} be divided into n disjoint subsets of equal size to, 
{!,..., N} — [J2=i such that the coupling constants within each Av 
are < 0, but > between Ai, and A^ for v ^ [i. Moreover, the partial row 
sums are assumed to be constant: 

T J.- 1^" '^"^-^'^ (29) 
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A spin system satisfying the assumptions of this section will be called n- 
homogeneous, see PI2||S|- Note that this notion is incommensurable to n- 
cyclicity defined in the previous section. However, certain rings, the trian- 
gular lattice, the kagome lattice, and the icosidodecahedron satisfy both 
definitions. A necessary condition for nearest neighbor Heisenberg systems 
to be n- homogeneous is that the number of nearest neighbors, which is as- 
sumed to be constant, is divisible by (n — 1). Actually, spin rings of even N 
are 2-homogeneous, rings of odd N are 3-homogeneous if N is divisible by 
3. n- homogeneous Heisenberg rings do not exist for n > 3 because they do 
not fulfill the homogeneity condition (|29|l . 

We recall that 1 = -^{^, 1, ■ ■ • 1 1) is an eigenvector of J with eigenvalue 
j. Due to n- homogeneity there are, after a suitable permutation of the spin 
sites, further eigenvectors of the form 



i^'') = (m : 1,771 



p^77^ 



„2k 



(n-l)fc\ 



1 



(30) 



where (777, 



.) denotes the Tn-fold repetition of the entry x, and p 



The corresponding eigenvalues are — j™ -\- J2p=i P' 



pk 



hence they coalesce into one (77 — l)-fold degenerate eigenvalue. By 
applying the theorem of Gersgorin (c.f. IT], 7.2) this eigenvalue is shown 
to be the smallest one jmin- _ 

Next we construct a coupling matrix JJ with the same eigenspaces as J 
but only three different eigenvalues. It has the block structure 

[ AC C ...\ 

C AC .. 

CCA... . (31) 



V 



where A and C are 777 x 777-matrices of the form 



A 



—a (3 ~a 
—a —a (3 



V 





/7 7 7 . . .\ 




777... 




777... 


J 





The three eigenvalues of JJ are 

J = /3 — [m — l)a + {N — 777)7, 
Jmin = /3 - (777 - l)a - 7717, 
j2 = a + ^ 

with degeneracies 1,77 — 1 and N — 11, resp. By choosing 

a = ;^("v?2 - i-{n- l)jmin), 

/? = ^((^ - n)j2 +J + {n- l)jmi„), 



7 = 



N 



(32) 



(33) 
(34) 
(35) 

(36) 
(37) 
(38) 
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one obtains 

J = j, Jmin = jmin, h = 32 ■ (39) 

j2 is the remaining smallest eigenvalue of J' after eliminating [n— l)-times 
jmin from the set of eigenvalues. Thus it can happen that j2 — jmin if jmin 
is more than (n — l)-fold degenerate. 

Let us write = for any subset A <Z {1, . . . , N}. We conclude 

H>H = -a{^ + (a + /3)iVs(s + 1) + 7 J^S^ - ^ S^,^^ (40) 

N N 

> jS{S + 1) - (a + 7)n— s(— s + 1) + (a + P)Ns{s + 1) . (41) 
n n 



3.3 Analytical expression for the lower bound 

Hence we obtain for the lower bound 

E > l^^^SiS + 1) + iVj„,ins(s + 1) + (TV - n){j2 ~ Un)s . (42) 

Since j2 — jmin > the bound (|42|l is the better, the smaller n is. This is 
in contrast to the upper bound considered in the previous section, which is 
improved for large odd n. 



4 Examples 

In the following examples we calculate the energy eigenvalues by numerical 
methods as well as lower and upper bounds. All examples are Heisenberg 
spin systems where the total spin 5 is a good quantum number. It turns 
out that S t—^ i?min(5') is always a monotonically increasing function, hence 
we need not to distinguish between -E'min('S') and E^in{M)- 

In order to judge the quality of the bounds we provide the deviation of 
the best upper and lower bound from the exact ground state energy in rela- 
tion to the energy difference between antiferromagnetic and ferromagnetic 
ground state, i. e. 

_ |^bomid,0 — Eo\ 

" ~ E{Ns) -Eo ' ^ ' 

The first example we would like to consider is a Heisenberg spin ring 
with iV = 6 and s = 1/2 as well as s = 5/2. Figure [21 shows the numeri- 
cally determined energy eigenvalues (dashes) as a function of total spin S. 
The solid curves display the bounds for the minimal energies considering 
2-cyclicity and 2-homogeneity. 

As a second example we take a frustrated Heisenberg ring with = 9 
and s ~ 1/2 as well as s — 3/2. The results are presented in Fig. 13 The 
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Fig. 2 Upper and lower bounds of -Emin(5') for Heisenberg spin rings with A'^ = 6 
and s = 1/2 (top) as well as s = 5/2 (bottom). The solid curves display the 
bounds for the minimal energies considering 2-cyclicity (s = 1/2 : e = 0.16; s = 
5/2 : e = 0.05) and 2-homogeneity (s = 1/2 : e = 0.28; s = 5/2 : e = 0.07). 



solid curves display the upper bounds for the minimal energies consider- 
ing 9-cyclicity, the dashed curves do the same for 3-homogeneity. Without 
using the concept of n-homogeneity the lower bounds are much poorer for 
frustrated systems [H]. 

Another example, an icosidodecahedral Heisenberg spin system, is re- 
lated to magnetic molecules, which can be synthesized in such structures. 
One species is given by {Mo72Fe3o}, a molecule where 30 Fe^+ paramagnetic 
ions (s = 5/2) occupy the sites of a perfect icosidodecahedron and in- 
teract via isotropic nearest-neighbor antiferromagnetic Heisenberg exchange 

. Not much is known about the spectrum of such giant structures since 
the Hilbert space assumes a very large dimension of 6'^° ~ 10^'^. So far only 
DMRG calculations could approximate the minimal energies |14|. 
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N=9, s=1/2 




12 3 4 
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Fig. 3 Upper and lower bounds of Emin{S) for Heisenberg spin rings with A'^ = 9 
and s = 1/2 (top) as well as s — 3/2 (bottom). The solid curves display the 
upper bounds for the minimal energies considering 9-cyclicity (s = 1/2 : e = 
0.20; s — 3/2 : e — 0.09), the dashed curves do the same for 3-homogeneity 
(s = 1/2 : e = 0.26; s = 3/2 : e = 0.10). 



Figure 01 shows as dashes on the l.h.s. the minimal energies for s = 1/2 
which are determined numerically by J. Richter with a Lanczos method |15l 
[T^ and on the r.h.s. the minimal DMRG energies The icosidodeca- 
hedral Heisenberg spin system is 3-cyclic as well as 3-homogeneous. The 
corresponding bounds are displayed by solid curves. Especially the upper 
bound for the case of s = 5/2 is very close to the "true" (DMRG) minimal 
energies and thus could be used to justify approximations of the low-lying 
spectrum as used in Ref. ^7|. The lower bounds are worse than expected, 
but this behavior is explained by the 10-fold degeneracy of jmin, therefore 
j2 — jmin, and the last term in ()42(l yields zero, unfortunately. 
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s 
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Fig. 4 Upper and lower bounds of _Bniin(<S') for Heisenberg spin systems with 
icosidodecahedral structure, i. e. A'^ = 30 and s — 1/2 (top) as well as s = 5/2 
(bottom). The solid curves display the bounds for the minimal energies considering 
3-cyclicity (s = 1/2 : e — 0.18; s — 5/2 : e = 0.03) and 3-homogeneity (s = 1/2 : 
e = 0.36; s = 5/2 : e = 0.10). 



The last example discusses the triangular spin lattice which is one of the 
frustrated two-dimensional spin systems. The triangular spin lattice is 3- 
homogeneous and 3-cyclic, if the periodic boundary conditions are suitably 
chosen. Figure |31 displays the energy levels for = 12 and s = 1/2 (l.h.s.) 
as well as s = 1 (r.h.s.). The bounds of i?min(<5') are given by solid curves. 
In both cases the upper bound is very close to the exact minimal energies. 

For the thermodynamic limit A'^ oo of the triangular lattice with 6 = 1 
we rewrite the bounds by introducing a continuous spin variable Sc = S/N 
running from to s and using j = 6, jmin = — 3 with twofold degeneracy, and 
limjv^oo j2 — —3. After dividing by N the resulting bounds are separated 
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Fig. 5 Upper and lower bounds of EminiS) for the triangular spin lattice with 
TV = 12 and s = 1/2 (top) as well as s = 1 (bottom). The solid curves display the 
bounds for the minimal energies considering 3-cyclicity (s = 1/2 : e = 0.09; s — 
1 : e = 0.05) and 3- homogeneity (s = 1/2 : e = 0.34; s = 1 : e = 0.19). 



only by 3s: 



- -3s< lim ^"""^^^ < - 



(44) 



A Proof of Eq. 

For easy readability we repeat Equation 
(f2M\A\[2M 



5\r\ f^^2_ 2sai2Ns-a) 



N 



2Ns-l 



(45) 
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Since 1?m is invariant w. r. t. arbitrary permutations of spin sites it 

suffices to choose A = s^^ and to multiply the result for (I^mI^I^^m) 
by (5 1 r|. We note that 

[A5-] =-(sf^S2 +sr4'^), (46) 

and 

[[A,S-],S-]=2s^s^, (47) 
but higher commutators vanish. Hence 

[A, {S-r] = a{S-r-' [A, S-] + Q iS-r-' [[A S-] , S-] . (48) 

Further we define A(a, A:) by 

(5+)"(5-)"(5-)'=|/?) = A(a,fc)(5-)'=|/2) . (49) 
Using S~^S~ = S"^ — S'^^^S^^^ — 1) one derives the recursion relation 

A(a +l,k) = {2Ns -a- k){a + k + l)A(a, A:) . (50) 
Together with A(0, A;) = 1 it can be solved and yields 

Xia k) - (2iV^-fc)! ja + ky. .5^. 
^^'''''>- {2Ns-a-k)\ k\ ■ 



Obviously, 



CM-A(a,0) - ^2Ns)\a\ ' ^^^^ 



hence 



C^\[a - 1, 1) - ^27Vs)! a! {2Ns - a)! 1! " 27V^ ' ^^^^ 



and 



GMA(a 2, 2) - ^2^^^, ^, ^2iVs - a)! 2! " 4Ars(2iVs - 1) ' ^^^> 
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Now we are prepared to calculate (i7A/|Z\|i7M): 

{[2M\A\nM) = cI,{Qm\{s-ta + [A, {s-y] \[2) (55) 

= + aClt{UM\{S-T-^ [A.S-] \n) (56) 
+ {^^Cl{QM\{S-r-^[[A,S-],S-]\Q) 

= s2 + aCl[{{s+f-^{s-T-^s- n\ [z\, S-] \n) (57) 
+ {^^cU{s^T-\s-r-^{s-fQ\ [[AS-] ,s-] \n) 

= + aCl[\{a - 1, 1){S-Q\ [A, S-] \n) (58) 
\Clj\(a-2,2){{S-fQ\ [[A,S-],S-] \n) 



-^+2^^-^-^)H2j 4iV.(2^.-l) ^^^' ^''^ 

AT I OAr„ II ' 



N y 2Ns - 1 

In line (|^ we used (@Hl- is obtained by means of (|^ . (|^ and the 

identities 

{s-n\ [A, S-] \n) ^ {s-n\ - (sf ^s" + 574^^)1^2) (61) 

= ~2(s2"i7|sf^S2"|J7> = -2s(s(s+ 1) - s{s - 1)) 
= -As' , 



and 



{{s-fn\ [[A,s-] ,s~] \n) = {{s-)'n\2s^s^\n) (62) 

= 4:{s^s^n\s'^s2\n) 

= 4(s(s + l)-s(s-l))2 



This completes the proof. 

Acknowledgements We thank Johannes Richter for fruitful discussions. 
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